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Abstract 

We study the stability of a contracting silent universe, which is a spacetime with 
irrotational dust and vanishing magnetic part of the Weyl tensor, = 0. Two 
general relativistic backgrounds are analyzed; one is an attractor of silent universes, 
i.e., a locally Kasner spacetime, and the other is a particular class of inhomogeneous 
Szekeres solutions. In both cases their stabilities against perturbations with non- 
zero magnetic part depend on a contraction configuration; a spindle-like collapse is 
unstable, while a pancake-like collapse is still stable. We also find that a similar in- 
stability exists in spindle collapse for a Newtonian case. We conclude that H a i, = 
is not a generic ansatz even in general relativistic dust collapse. 
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1 Introduction 



The recent observations of the large scale structure in the universe show that quite non- 
linear wall-like structures or filament-like structures seem to exist. Such structures may 
give us important information about the structure formation process in the universe. It 
may be necessary to study nonlinear dynamics in order to know what kind of structure 
is formed on large scales. 

In Newtonian gravity, the Zel'dovich approximation (Zel'dovich 1970), which describes 
well the evolution of inhomogeneities even in the weakly nonlinear regime, predicts a 
pancake collapse. The analysis of a homogeneous ellipsoid in an expanding universe also 
supports such a pancake collapse (White & Silk 1979). How about for a very large scale 
such that general relativistic effects become important? 

In general relativity, however, the analysis of nonlinear dynamics is quite difficult 
because the Einstein equations are highly nonlinear partial differential equations. It has 
been restricted to very simple cases such as the spherically symmetric Tolman-Bondi 
model (Tolman 1934; Bondi 1947). Recently the nonlinear dynamics of irrotational dust 
with vanishing magnetic part of the Weyl tensor (H a b = 0) has been examined (Barnes & 
Rowlingson 1989; Matarrese et al. 1993; Croudace et al. 1994; Bertschinger & Jain 1994, 
Bruni et al. 1995a, b). Under such a condition, the system is described by nonlinear 
ordinary differential equations, which makes us able to analyze more generic situations. 
Such a system has no degree of freedom of gravitational waves, no information can transfer 
between the neighborhoods, and each fluid element evolves independently. Therefore such 
models were called silent universes by Matarrese et al. (1994). Bruni et al. (1995a) found 
attractor solutions in that system and showed that a generic gravitational contraction 
leads to a triaxial spindle-like configuration. 

However the assumption of H ab = is controversial. It may not be necessarily justified 
for generic cases. In fact, Matarrese et al. (1993) and Kojima (1994) pointed out, using 
a perturbative approach in an expanding Friedmann-Robertson- Walker (FRW) universe, 
that even if H ab = in a linear regime, H ab will be generated during the nonlinear 
evolution of inhomogeneities and would affect the dynamics of the universe. This leads to 
another question, i.e., "Is the attractor in the silent universe still some attractor in more 
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generic spacetimes?" 

In this paper, we study perturbations of silent universes, taking into account the effect 
of non- vanishing H ab . We shall call such a perturbed spacetime a quiet universe. Then 
the question above becomes "Will a quiet universe be silent ?" Analyzing stability of 
the attractor solution of silent universes (Kasner-type spacetime) and some particular 
class of Szekeres solutions (Szekeres 1975), we discuss whether or not the silent universe 
is generic in general relativity and what kind of configuration is more probable in the 
nonlinear regime of structure formation. 

This paper is organized as follows. In §2, we write down the basic equations for the 
dynamics of irrotational dust with the introduction of dimensionless variables and we 
summarize the silent universes. In §3, the stability of silent universes against linear per- 
turbations with non-vanishing H ab is analyzed. We consider two background spacetimes 
as unperturbed spacetimes: the Kasner-type universe and some Szekeres universes and 
show the dependence of stability on contraction configuration. We also discuss the sta- 
bility of the Zel'dovich solution in Newtonian theory in §4. Our conclusion follows in 
§5. 

We take c = 87rG = 1, and the signature of (— , +, +, +). Throughout this paper, Latin 
indices are used for coordinate components; a, b, c, ■ ■ ■ run from to 3 and i, j, k, ■ ■ ■ are for 
the three spatial indices, while Greek indices are used for tetrad components; a,j3, 7, • • • 
run from to 3 and X,fj,,v,--- are for spatial triad indices. The brackets () and [] denote 
symmetrization and anti-symmetrization for all indices. 

2 Dynamics of Irrotational Dust 

2.1 Basic Equations for Irrotational Dust 

The irrotational dust fluid system with four-velocity u a is described by the following 
variables: the mass density p, the expansion scalar 6, the shear tensor a ab and electric 
and magnetic parts of the Weyl tensor C abc d] E ab and H ab . Those are defined as 

9 = u c ]C , (2.1) 
(Tab = U( a;b ) - -h ab 9, (2.2) 
E ab = C acbd u c u d , (2.3) 
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H ah = *C acbd u c u d , (2.4) 

where h ab = g ab + u a u b is a projection tensor, and *C abcd = \^abefC e ^ cd is the dual of the 
Weyl tensor. The basic equations (Ellis 1971) are 

+ U 2 + 2a 2 + \p = 0, (2.5) 
p + 6p = 0, (2.6) 

2 2 2 

c^b + o" ac o" C 6 + -^Qvab - ^h ab (T + E ab = 0, (2.7) 
h ac h bd E cd + h c(a e b)def u d H ce .j + h ab a cd E cd + 6E ab - 3E c {a a b)c 

= ~\p° a \ ( 2 -8) 
h ac h bd H cd - h c{a t b)def u d E ce . f + h ab a cd H cd + - 3# C ( V 6 ) C 

= 0, (2.9) 

V-< fe = 0, (2.10) 

H ab = -h ac h bd cr e ic ]f e d ^ ef u g , (2.11) 

hTh^E^ - e abcd u b a e c H de = h ab p, b , (2.12) 

h ac h bd H cb , d + e abcd u b a\E de = 0, (2.13) 

where, an overdot represents a directional derivative along the fluid worldlines, e.g., T a ... b = 
T a ... b ; C u c . 

We introduce the orthonormal tetrad {e^} by 

e (0 a ) = u a = (1, 0, 0, 0), e ( ; } = (0, e^), (2.14) 

with 

e^e? = 6;. (2.15) 
The spatial triad vectors {e^} are chosen to be Fermi-propagated along u a , i.e., 

e ( %u b = 0. (2.16) 

The metric is described by those tetrad vectors as, 

ds 2 = g ab dx a dx b 



-dr + e^e^dx'dx 3 . (2.17) 



To write down the basic equations in terms of the tetrad components, we use the following 
relation for the covariant derivatives: 

= — T 5 ia Ts/3 — r 5 7/3 T a 5, (2-18) 

where T a/ 3 y = e^e^a^e 6 ^ is the Ricci rotation coefficient and d a represents the direc- 
tional derivative with respect to e( Q ). The governing equations for the tetrad components 
then become the same ones for the coordinate components in Ellis (1971). 

We also need the equations for the Ricci rotation coefficients r°jg or the structure 
constants 7°^ = r°jg — T a yf3 to close our system. Those equations are obtained from the 
Jacobi identities, i.e., 

%T a ^ + 7[^TV = a ( 2 - 19 ) 
From our choice of the tetrad, some of the components of r°jg 7 are known as 

r° 00 = o, n = o, r° 0u = o, r" „ = o (2.20) 

and 

7^o = n o = ^ + ^M- (2.21) 
Next, we define dimensionless variables, Q, E^, 7"^,, C/^p and C M ^ P as, 

P = g^^ 2 ' ° iiv — E^jA i?^ = S^O 2 , H^ u = Ti.^9 2 , 

r^^p = Q^vp^i Ifxvp = ^tavp^i (2.22) 



and 



E 2 = is pCT E^. (2.23) 



Using those variables, an attractor solution in silent universes becomes a fixed point in a 
"phase" space. Now the basic equations for irrotational dust fluid are 

9' = i(2 + fi + 12E 2 )fl, (2.24) 
ft' = i(l -n- 12E 2 )fi, (2.25) 
E M ; = ^(2-Q-12E 2 )E^ + E w E^-^E 2 + £:^, (2.26) 
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+ 2^(^V CT + e MP^/ ;CT ) + ^{e vpa U p + e w «H P X (2.27) 
?C = ^(l-fi-12S 2 )^ + 5^ p(T S^-^(^ p S^ + S w ^) 



2^ {tvpa^if' + £[ip<j£ u P ' ) ^2^' ^ tu P a ^ P p e P*pv£ fP v)i (2.28) 



^(E^-^) + S^ = 0, (2.29) 
= "^(^S/ + ^/-) - ^(^/ + ^S/-), (2.30) 
V + \** B v> = e»»eV"H°» + ^ + §W (2.31) 

+ 2 f' u n»v = -e^e xra* , (2.32) 

where prime (') denotes the derivative with respect to r = — / 6dt and e pup is the Levi- 
Civita tensor, r ~ ln(a proper volume) -1 plays the role of time in the contracting phase. 
The Jacobi identities are 

d [s c^ + l -c a [5l d p] e + e c\ Sl c% = 0, (2.33) 

which consist of 12 independent equations: 9 dynamical and 3 constraint equations. Hence 
we have now 26 evolution equations and 17 constraint equations. 

2.2 Silent Universes 

The silent universe is proposed by assuming that (i) the fluid is collisionless dust (p = 0), 
(ii) it has no vorticity {uo^ = 0), and (iii) the magnetic part of the Weyl tensor vanishes 
(Hob = 0). Such an ansatz seems to be plausible because gravitational waves may have 
nothing to do with the structure formation process^. 

a a b and E ab are diagonalized under this condition as shown by Barnes & Rowlingson 
(1989). The two independent components are diagonal: an, <j 2 2 and En , E 22 - The 
metric is written as 

3 

ds 2 = -dt 2 + Y,ll{x 1 t){dx t f, (2.34) 



1 It turns out, however, that a contribution of the magnetic part of the Weyl tensor may not also be 
negligible even in the Newtonian limit (Bertinger and Hamilton 1994). 
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and, 



(2.35) 



where the Z^'s represent the scale factors in the i-direction and an is the i-i tetrad com- 
ponent of the shear tensor. 

For the silent universes, the evolution equations are reduced to the following closed 
set of ordinary differential equations. 





n 

E_ 
E_ 

£ 



^(2 + fi + 12E 2 ) 9, 



12E 2 n 



6 
1 

3 

- (2 - o - 12E-) e 
1 



s 4 



QK 2-Q- 12E 2 J E_ + + £-> 

- (l - n - 12E 2 ) £+ + v^S^ - s_£_ 

- (l - O - 12E 2 ) £_ - y/3(E + £_ + E_£ 4 



o 
6 



(2.36) 
(2.37) 
(2.38) 
(2.39) 

(2.40) 
(2.41) 



whereQ 



S 4 



V3 
2 

^3 
2 



(En + E22), 
(£11 +£22), 



2n 



:(^ii-^ 



J 22J 



22 y 



and 



Ei + E^ 



(2.42) 



(2.43) 



In the present case, the equations for Q a /3-y, or C a /? 7 are decoupled. From the condition 
of = 0, we also find a new constraint equation. However, those are guaranteed to be 
satisfied from the dynamical and other constraint equations (Lesame et al. 1994). 

Because of the absence of spatial derivatives in the evolution equations, each fluid 
element evolves without influence of the neighborhood. Hence it is called a silent uni- 
verse. It includes not only homogeneous Bianchi I spacetimes, but also inhomogeneous 
spacetimes such as spherically symmetric Tolman-Bondi spacetimes or Szekeres solutions. 



2 This choice of ± variables is different from the definition by Bruni et al (1995a), but is the same as 
that in usual discussion on Bianchi models. 
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We see the evolution equations for 5 dynamical variables {O, E±, S±} are closed by 
themselves. Analyzing them, Bruni et al. (1995a) find attractor solutions in 5-dimensional 
"phase" space of {Q, £±, S±\. 6 is obtained from eq. Q2.36|) after solving those 5 variables. 

The attractors in contracting spacetimes (6 < 0) locally correspond to a family of 
vacuum Kasner spacetimes. The set of attractors is a circle in a 2-dimensional £_} 
"phase" space (Fig. 1). Only three points on the circle of attractors correspond to a 
pancake collapse, i.e., two of the U stay finite but the other one vanishes. The rest of the 
points of the attractors correspond to spindle collapse, i.e., two of Zj's vanish but the other 
one diverges. Bruni et al. (1995a) then concluded that contracting regions generally tend 
to be spindle-like configurations in silent universes. 

3 Quiet Universes: Magnetic Perturbations of Silent Universes 

In order to generalize the analysis by Bruni et al. (1995a), we consider linear perturbations 
with H a b of a silent universe. First we study the case in which the background solution 



is an attractor of contracting silent universes, i.e., the Kasner-type spacetimes in § pTl 
Secondly, as a background solution, we adopt a planar inhomogeneous Szekeres solution 
which approaches a particular attractor solution of silent universes (§ |3.2| ). We denote the 



background quantities of the silent universe by a bar ( ) and the perturbed ones with tilde 



and 



e = e(i+#), n = n + n, s^ = s^ + s^, 

£flV £[1V + ^flUJ H/AU 'H/lVi 

Gfivp Gfivp + Gpvpi C^pp Cpiup + Cp,yp, (3.1) 



S 2 = ^.XT. (3.2) 



Note that = because the background spacetime is a silent universe. 

The background spacetime we analyze here has at least one homogeneous direction, 
which we choose to be x (or 1). We then analyze only a single plane- wave mode of 
perturbations with a comoving wave number k = (k, 0, 0). The spatial gradient d\ is 
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replaced as 



(3.3) 



In this case, the perturbation equations, which consist of 26 dynamical and 17 con- 
straint equations, are divided into the following 4 groups, where we have introduced {E±}, 
{E±}, {£±}, {£±}, and {H±} in the same way as Eq. (|2.42|) . 

(1) Group 1 (0, fl, E±, 4, 7i 23 , C 33 i, C 221 ): 
(1-1) 9 dynamical equations: 



7? 



ft 



23 



-(2 + + 12S 2 ) 5 + ^fi + 4(E + E+ + E_E_), (3.4) 
6 6 

hi -si- i2s 2 )^^ + -(1 - 2n - 12s 2 ) n 

3 3 

-80(S + S + + S_S_), (3.5) 
is + (2 - O - 12S 2 ) - -^(E 2 , - E 2 ) + 4 £ - 

+-(2 - tt - 12E 2 )E+ - 4E + (E+E + + E_E_) - 4=(E + E + - E_E_) 
6 ' V3 

-4, (3.6) 



(2 - n - 12s 2 ) + 4=e+e_ + £_ 

6 v 



>/3 



6 



+-(2 - Q — 12E 2 )E_ - 4E_(E + E + + E_E_) + 4=( S +K 
6 y3 

+4, 



E_E + ) 



(3.7) 



- fl - n - 12E 2 ) £ + + V3(Z + £ + - E_£_) + is+fi 
3 v 7 6 

+i(E+-24)n 

6 

+hl — Ct — 12E 2 )4 + >/3(S+4 - S_4 

o 



1? 



6 

v/3, 



-(8E+ - v / 3)4S+ " (8E_4 + \/3£_)E_ + -^-0^23, 
- (l - Cl - 12E 2 ) £_ - \/3(E + £_ + E_4) + 5 
+i(E_ -2£_)fi 

+-(1 - fi - 12E 2 )£_ - v^fS+4 + e_4) + 
3 6 

-(8E+ + v / 3)^-S + - (8E_£_ + V34)E_ - 



(3.f 



(3.9) 



i(l - n - 12E 2 ) + ^(E + + v^S^) 



^23 + (V3S+ - £-)di# 
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+i<9 1 (v / 3£+ - £-) - Z-Q22X + \(VZ£ + + £-)G 3 si, 



r - 

u 331 — 



a 



221 



In + 2S 2 -i=(S + + v / 3S_) 
-±(l-2y/3X + )d 1 #, 

In + 2s 2 - --*=(e + + v^s^) 

v3 
1 + _L (S+ -V3E.) 



^331 + yg^ 1 ^- 



(1-2) 3 constraint equations: 



[(£+ + v 7 ^) - -=] + 9i(E + + v^X-) + v^v^ + E_) C. 
V3 



331 



+ 2v / 3S_C 221 = 

= - XL) + 9 1 (v / 3S+ - XL) + (V/3S+ + XL) C ; 

— 2E_C 22 i], 



331 



+ 2[(£+ + V^5_) 



3^ 



v / 3(v / 3S + - E_)7i 23 + ^jf^. 



(2) Group 2 (E 23 , £ 23 , H±, C 123 , C 231 , C312) : 
(2-1) 7 dynamical equations: 



E 23 = -[2-n-12E 2 -2v / 3(E + + v / 3E_)]E 23 + ^ 23 , 



'23 



n 4 



n_ 



'123 



6 



-(1 - Cl - 12E 2 ) + — (E+ + v^E^) 



'23 



+ hn + 3V3(£+ + Vs£-)]t 23 - ^d 1 (Vsn+ - it-), 

1 _ ~ /3 

-(1 - ft - 12E 2 )7i: + + >/3(E+W+ - E_7i_) - ^-9x^3 

-^[(>/35 + + £-)Qi2z + (v^ - B-W312], 
i(l - ft - 12E 2 )7lL - a/3(£+W- + E_7i+) + ^^ 23 
- ^(VZ£+ + £-)Gi23 + \{^>£+ ~ £-%i2 ~ 25-&31, 



l0 + 2E 2 + -|(E + + v / 3E_) 



-123, 



c. 



231 



-n + 2s 2 + 4=(s+ - 73S_ 

6 y/3 K + 



^231 — <9iE 23 , 



^312 — 

(2-2) 4 constraint equations: 



^+ = ^-[<9l£ 2 3 + ^ 12 3 ~~ ^- ^231 + V^E^- C312], 

= -[^E.g - (V/3S+ + 2E_) C 123 + (V/3S+ - 2E_) C 23 i 
+ S-C312], 

d 1 (H + + v^ft-) = >/3(>/35 + - £-)E 23 - V3(Vsi: + - s_ 
diC 123 = 



(3) Group 3 (E 31 , £ 31 , H 12 , C 223 , C 113 ): 
(3-1) 5 dynamical equations: 



J 31 



6 



5. 



31 



^12 = 



--223 



'113 



[2 - n - 12E 2 - 2^(2+ - v^E-)]^! + £ 3 i, 
i(l-n-12£V^(£+-v / 3£-) £31 

+ I[Q + 3V3(£ + - V3S-)]t 31 - U W 12 , 
o 2 

i(i - n - 12s 2 - 3 V / 3s+)7i: 12 - ^£31 

+ \[(V3£+ + £-)&i 3 - (Vs£ + - £-)Gt23], 
-(^ + 2E 2 + ^E + )c 223 , 

- Qft + 2E 2 + -^=E + ^ C113 + 9iE 3 i, 



(3-2) 5 constraint equations: 



- (V/3S+ + S_) C113 - (v/3S + - E_)C 223 = 0, 
H 12 = -i[9xS 31 + (v/3S + - S_) C 223 - (v/3S + + E_) C 113 ] 
9i5 3 i - (VS£+ - £-) C223 - (VS£+ + £_) C113 = 2E_H 12 , 
d 1 H 12 = -(V3£ + + £_)E 31 + (V/3S+ + E_)£ 3 i, 
3iC 223 = 
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(4) Group 4 (E 12 , £ 12 , H u , C 112 , C 332 ): 
(4-1) 5 dynamical equations: 



J 12 



g ' - 
°12 ~ 



-(2 - ft - 12E 2 + 4v / 3E + )E 12 + e 12 , 
o 

-{1-Cl- 12E 2 - 3v / 3S+)^i2 - (V3S+ - -ft)E 
3 6 



12 



31 



I(l-ft-12E 2 ) + ^(E + -v / 3E_) 



112 



- -(y/3S+ - £-)Q 3 32 - £-Gll2, 

-Cl + 2E 2 - 4=(£+ - v^S- 



332 



ft + 2E 2 - 4=fS. 



Cll2 + 9iEi2, 
^332? 



(3.37) 
(3.38) 

(3.39) 
(3.40) 

(3.41) 



(4-2) 5 constraint equations 



<9iE 12 + (V3Z + - E_) C 332 - 2E_C 



112 



0, 



W 3 i = 2^12 " 2S_C 112 - (v/3E + - S_) C332], 
di£n + (VS£ + - £_) C332 - 2£_C 112 = + £-)W 3 i, 

d\H 3 \ = 2£^_Si2 — 2E_£^i2, 
= 0. 



(3.42) 
(3.43) 
(3.44) 
(3.45) 
(3.46) 



Since each group is decoupled, 6 (in Group 1) and 3 (in Group 2) dynamical variables 
remain free at the initial time. Although we have fixed our tetrad system in the back- 
ground spacetime, we still have some freedoms for their choice in perturbed spacetime. 
The time is chosen as the proper time of a dust particle, while the spatial coordinates have 
been left to be free in perturbations. As for the physical dynamical degrees of freedom 
of perturbations, we must have scalar and tensor modes, which are coupled except for 
the FRW spacetime background. Notice that we do not have a vector mode because of 
our irrotational dust. Since our perturbations are inhomogeneous only in the 1-direction, 
the tensor modes are decoupled into two parts; one is coupled to scalar perturbations 
(called even parity mode) and the other is decoupled (called odd parity mode). The per- 
turbations of Group 1 then consist of scalar and even parity tensor modes, which have 4 
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dynamical degrees of freedom. The other 2 variables are left still to be free because of the 
freedom of coordinate transformation in 2- and 3- directions. The perturbations of Group 
2 consist of 2 physically dynamical freedoms (odd parity tensor modes) and 1 freedom of 
coordinate transformation in 1-direction. There is no degree of freedom in Groups 3 and 
4, i.e., no physical perturbations. 

Initially, we will specify Cl, £+,£+ and 7i 23 in Group 1 (corresponding to 4 physical 
modes). We set C331 = C221 = 0, using remaining gauge freedoms, and determine other 
variables from the constraint equations. Then we analyze stability against the perturba- 
tions with 4 physically independent initial variables, (fl, E + , £ + , H23) = (1 + 2,0,0,0), 
(0, 1 + i, 0, 0), (0, 0, 1 + i, 0), and (0, 0, 0, 1 + i). In Group 2, we determine the initial data 
from 2 physical variables in the same way as Group 1. 

For the first case of our analysis (§3.1), the choice of x direction and a single wave 
mode does not cause a loss of generality, because x is not a special direction in the present 
problem and any linear perturbations can be expanded by Fourier series. As for the second 
case (§3.2), the perturbations in another homogeneous direction (y) will have the same 
properties as those in the x direction. However, for the perturbations in the direction 
of inhomogeneity (z), we may not be able to define invariant perturbations because the 
background itself is inhomogeneous. Therefore we have not analyzed it here. 

3.1 Perturbations of the Attractor Solutions 

The attractor solutions of silent universes are homogeneous vacuum Kasner type space- 
times, which are characterized by a constant value of 77, (0 < 77 < 2tt), as 




(3.47) 



(3.48) 



£+ = — — (2 COST/ - 1)(1 + COS77), 




(3.49) 



£- = (2 cosr? — 1) sin 



(3.50) 



9 



and 



e = 



i + 0.(t-**)' 



(3.51) 
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where is the expansion at an initial time t — t*(T — r*). r is given as 

t-t* = -ln[l + < ,(t-t»)]. 



(3.52) 



We have analyzed the behaviors of perturbations for various background values (£+) 
and for several scales of the perturbations (A: -1 ). 

The typical behavior of the perturbations is shown in Fig. 2, in which we set the 



the other 3 independent initial data in Group 1, and find similar results. 

The evolution of perturbations is quite different depending on whether the background 



and (—-^, 0). In the case of a spindle-like collapse, the magnetic part of the Weyl tensor 
gets larger and larger, and eventually diverges. The "silence" is broken. Other variables, 
e.g., some components of shear and the electric part of Weyl tensor also diverge. This 
asymptotic behavior is qualitatively independent of the wave number k and background 
parameter r/. On the other hand, in the case of a pancake-like collapse, the perturbations 
do not grow for any wave number. In pancake collapse, we have also analyzed another 
two independent perturbation modes, i.e., Group 2 modes and confirmed that it is stable. 

As a result, we conclude that against generic perturbations with H^, a pancake-like 
attractor is still stable, but a spindle-like attractor becomes unstable. 

3.2 Perturbations of the Szekeres Solutions 

In the above analysis, since the background spacetime is a homogeneous attractor solution, 
our result might not be generic. It is possible that the evolution of the perturbations 
in an inhomogeneous background is different from that in a homogeneous background. 
Therefore we consider the following particular class of Szekeres solutions (Szekeres 1975; 
Goode, S. W., & Wainwright, J., 1982) as an inhomogeneous background silent universe: 



where (5{z) represents an inhomogeneity in the z-direction. When /3(z)—0, the spacetime 
is just a flat Friedmann-Robertson- Walker (FRW) universe. 



initial values of the perturbations to (Q, S + , £ + , H23) = (1 +h 0,0,0). We also analyze 



collapse is spindle-like or pancake-like. Three pancake points are (£+,£_) = (-4=,±| 





(3.53) 
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The dimensionless variables defined above for this solution are given as 

* - f~f, ( t (3-54) 

t(t-p(z)y 1 ; 

= 1-3E+, (3.55) 

£ + = -j=E + (l + V3E+), (3.56) 

£_ = s_ = 0, (3.57) 

where 

E + = r M . (3.58) 

t is given as 

T T - l n *1 Z gW} (3 59) 

There is no stationary point in the "phase" space (see Fig.l). In the case of contracting 
phase, when t decreases, even if the Szekeres universe is located initially near the isotropic 
FRW point (E+ = 0), such spacetimes fall into different attractors according to the sign 
of E + , i.e., when —l/y/3 < E + < 0, it falls into a pancake-like attractor (E + = — 1/\/3), 
while if < E + < l/y/3, then it approaches a spindle-like attractor (E + = l/\/3) (see 
Fig.l). 

Since in the z direction, the background Szekeres solution is inhomogeneous, we focus 
only on the perturbations in the x, y directions. Here we analyze only a single plane-wave 
mode with a wave number k = (k, 0, 0), as mentioned before. 

The results are shown in Figs. 3, in which we set the initial values of the perturba- 
tions (Cl, E + , £ + , H23) — (1 -M, 0,0,0). We also analyze the other 3 independent initial 
data in Group 1, and find similar results. For a spindle-like background (E + > 0), the 
perturbation of the magnetic part of the Weyl tensor grows infinitely and the "silence" is 
broken. The electric part and the shear also diverge. On the other hand, for a pancake-like 
background, the perturbations do not grow. The result depends only on the contraction 
configuration, i.e., the sign of E + . It does not depend on other factors such as the spa- 
tial scale of perturbations. In pancake collapse, we have again analyzed another two 
independent perturbation modes, i.e., Group 2 modes and confirmed that it is stable. 

Similarly to the Kasner-type background case, the present analysis shows that a 
pancake-like collapse is stable while a spindle-like collapse is unstable. One can also 
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find that an isotropic collapse (FRW) is unstable against perturbations with magnetic 
part H a b (what is called the tensor mode). 

Then some questions may arise, i.e., "Where will the destabilized solution approach in 
the "phase space" after a nonlinear evolution ?" or "How will the configuration change?" 
Our analysis shows that the perturbations with any wave length eventually approach the 
homogeneous one (k = 0). Therefore, we may expect that the spacetime always evolves 
locally into some Bianchi IX spacetime with complex oscillations near the singularity as 
shown in Belinskii et al (1970, 1982). 

4 Newtonian Case: Stability of Zel'dovich solutions 

Finally, in order to understand our results better, we analyze the stability of a silent 
universe in Newtonian gravity. The one-dimensional Zel'dovich solution in Newtonian 
gravity is very similar to the Szekeres solution in general relativity as pointed out by 
Kasai (1993), and then it may be regarded as a silent universe. Adopting the Zel'dovich 
solution as a background solution, we analyze its perturbations. 

In Newtonian gravity, we can describe the irrotational dust fluid dynamics similar 
to the case in general relativity except for the absence of the Maxwell-like equations of 
Eij and in general relativity (Ellis 1971). The system with irrotational dust fluid is 
described by a velocity field v l , mass density p, expansion scalar 9 = v 1 i; shear tensor 



and tidal force 



= 4>,ij - -5 i: j(j) ,l p 



where <fi is Newtonian gravitational potential. 

The evolution equations of p, 9, or have exactly the same expressions in general 
relativity, but E^ is determined only by its constraint equations without its evolution 
equation. does not appear. The equations for dimensionless variables, Q, S^-, S^, 
defined by the same equations as (|2.22|), are 



9' = l(2 + n + 12J: 2 )9, (4.1) 



n' = -(i-Q-i2i: 2 )Q, (4.2) 

3 
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= ^(2-fi-12S 2 )S ij + S^ i J -^S 2 + ^, (4.3) 



^l N :" • >: ,? , ; o. (4.4) 



e ilm (X jl;m + gO^i) + e jlm (X il ' m + -d> m H a ) = 0, (4.5) 
^J + \ B ^ = l^+l&Hi), (4.6) 
£ Um (£ji,m + - e 0, m £ 3l ) + e jlm (£ U ' m + - e e> m £ il ) = 0, (4.7) 

where £ 2 = E lm E lm /2 . 

The one- dimensional Zel'dovich solution is obtained by transforming from the Eulerian 
coordinates r = (r 1 ,r 2 ,r 3 ) to the Lagrangian coordinates q = (qi, q 2 , 93); 

ri = P /3 qi, (4.8) 
r 2 = t 2/3 g 2 , (4.9) 
r 3 = ^fe-UCt)^)], (4.10) 

where S(qs) is an arbitrary function of q^. Since we are interested in the contraction phase, 
we choose D(t) = t^ 1 as t decreases. This choice is different from the usual Zel'dovich 
solution, in which the growing mode of linear perturbation theory in an expanding universe 
is chosen (D(t) = t 2 ^ 3 in the case of Einstein-de Sitter model). Although the perturbation 
of the Zel'dovich solution was already studied by Bildhauer et al. (1992), we reanalyze it 
here because our choice of D(t) is different from theirs. 

Transmission of information will take place due to the Poisson equation, so non-locality 
exists in general as is well known in Newtonian gravity. However, the Zel'dovich solution 
can be called a silent universe, because once an initial condition is set, that is, S(q3) is 
given, each fluid element evolves independently of other neighbors. 

Therefore, the perturbation analysis can be done in the same way as in § |3.2j . In fact 
the background variables are given by the same solutions as eqs. ( |3.54| )~( |3.58| ) with 



replacement of f3(z) with S(qs). We analyze only perturbations in the x direction with a 
wave number k = (k, 0, 0). The perturbation equations are 

d = - (2 + Q+ 12S 2 ,) $ + -fi + 4S + S + , (4.11) 
6 v ' 6 
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- (l - O - 12 E+) tt # + - (l - 20 - 12 E+) 



8 

1 



6 



2 - - 2v / 3(l + 2v / 3S + )£+) E 



+ £+ + i[2-0 - 4v / 3S+(l + 373£_ 



6 

2-0 + 4v / 3S+(l - \/3E_ 



(>/3E+ - 2) ! + \/3(E+ + V3E-) t i = 0, 
3V3(£+ + v 7 ^-),! + 2(3V3£+ - 0) i = Oj, 
V3E+0, i + (v 7 ^ - = 0, 
£,i + - = 0, 

^12,1 = E13 I = E23 1 = £12,1 = ^13,1 = ^23,1 = 0. 



(4.12) 

(4.13) 
(4.14) 

(4.15) 
(4.16) 
(4.17) 
(4.18) 
(4.19) 



We have 6 variables with 4 constraint equations. Then the degree of freedom is two, 
which correspond to scalar perturbation modes. We find that the perturbations of Ejj 
and Sij behave very similarly to those in the case of the Szekeres background in general 
relativity. For a spindle-like collapse, the perturbations diverge, while for a pancake-like 
background, the perturbations decay (see Fig. 4, in which we set the initial values of the 
perturbations (Q, £ + ) = (1 + i, 0). We also analyze another set of independent initial data 
(Q,£+) = (0,1 + i), and find similar results). This result is consistent with the preference 
of pancake collapse in Newtonian theory. The results suggest that the instability of a 
spindle collapse shown in §|] is not due to an extreme situation near a singularity in 
general relativity. 

5 Conclusion 

We have studied the dynamics of a quiet universe, i.e., the perturbations of a silent uni- 
verse with non-zero H a b. The analysis of the silent universe by Bruni et al. (1995a) shows 
that a gravitational collapse prefers a spindle-like configuration. However, to consider a 
more generic case, we have analyzed perturbations with non-vanishing H a b for two gen- 
eral relativistic background solutions: (i) the attractors of the silent universe and (ii) a 
particular class of Szekeres solutions. In both cases, the behaviors of perturbations are 
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qualitatively the same. The stability clearly depends on the configuration of a contracting 
background spacetime, i.e., either a pancake-like or a spindle-like collapse. In the case of 
a spindle-like collapse, the perturbations diverge and the silence is broken. The space- 
time locally approaches the most generic homogeneous spacetime, that is, the Bianchi IX 
universe. On the other hand, the perturbations do not grow in the case of a pancake-like 
collapse. We conclude that even in general relativity, taking into account the effect of 
H a b, a spindle collapse is destabilized. 

This result agrees with the fact that a Newtonian collapse prefers a pancake-like config- 
uration. Our Newtonian analysis of the stability of the Zel'dovich solution also supports 
this result. To know the relation between general relativistic perturbations with non- 
vanishing magnetic part and the Newtonian ones, we study the behavior of the magnetic 
part in the relativistic equation ( j2.27| ) (see Hui 1995). We find that the terms with the 
magnetic part diverge for a spindle collapse but keep finite for a pancake collapse. This 
suggests that the instability for a spindle collapse in Newtonian theory may also be caused 
by the magnetic part which is responsible for transmission of gravitational information. 
We conclude that an ansatz of silence (H a b = 0) is no longer valid in generic gravitational 
contraction. 

We would like to thank P. Haines and Y. Kojima for useful discussions. This work 
was supported partially by the Grant-in-Aid for Scientific Research Fund of the Ministry 
of Education, Science and Culture (No. 06302021 and No. 06640412), and by the Waseda 
University Grant for Special Research Projects. 
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Figure Captions 



Fig. 1: A set of attractors of silent universes in the £_} plane; three dots represent 

pancake-like attractors. The rest points on the circle correspond to spindle-like attractors. The 
line of £_ = (— l/v3 < S + < l/v3) shows a particular class of Szekeres solutions adopted as 
the background solution in § |3.2| . 

Fig. 2: The evolution of perturbations for the attractor of silent universes. The initial values 
of the perturbations are chosen as (fi, E + , £ + , = (1 + i, 0, 0, 0). We set C221 = C331 = and 
other initial values are fixed by the constraint equations. Figs. 2(a) and 2(b) are in the case of 
pancake background (£+ = -1/V5). Figs 2(c), 2(d) and 2(e), 2(f) are in the case of spindle 
background (£ + = — 1/(4\/3) and l/\/3, respectively). The solid and dashed lines show the 
cases of the wave number k/lu0* = 10 and = 1, respectively. 

Fig. 3: Fig. 3: The evolution of perturbations for a particular class of Szekeres solutions. 
The initial values of the perturbations are chosen as (Cl, £+, £+, W23) = (1 + i, 0,0,0). We set 
C221 = C331 = and other initial values are fixed by the constraint equations. Figs. 3(a) and 
3(b) are in the case of pancake background (£+ = — 1/(4\/3))- Figs 3(c) and 3(d) are in the 
case of spindle background (S+ = l/(4\/3)). The solid and dashed lines show the cases of the 
wave number k/lud* = 10 and = 1, respectively. Figs. 3(e) and 3(f) are in the case of isotropic 
background (S + = 0). As for the isotropic (FRW) background, as was well known, the scalar 
perturbations are decoupled from the tensor modes. Hence, although our perturbation with 
the initial value (f2, £+, £+, W23) = (1 + i, 0,0,0) diverges and the FRW spacetime is unstable 
against this perturbation(Fig. 3(e)), W23 remains zero because it is a scalar mode. As for initial 
perturbations with H23 7^ 0, H23 diverges and then the FRW spacetime is unstable against this 
tensor perturbation as well(Fig. 3(f)). Fig. 3(e) has just a solid line because the perturbation 
equations for the scalar mode of isotropic dust universe do not involve the wave number. 

Fig. 4: The evolution of perturbations for the Zel'dovich solutions. The initial values of the 
perturbations are chosen as = (1 + i, 0). Other initial values are fixed by the constraint 

equations. Fig. 4(a) is in the case of a pancake background (£+ = — 1/(4\/3)), while Fig. 4(b) 
is in the case of a spindle background (£+ = l/(4v3)). Figs. 4 are drawn with just a solid line 
because the perturbation equations do not involve the wave number. 
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